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ABSTRACT 

We search for signatures of planets in 43 intensively monitored microlensing events that were observed 
between 1995 and 1999. Planets would be expected to cause a short duration (~ 1 day) deviation on the 
smooth, symmetric light curve produced by a single-lens. We find no such anomalies and infer that less 
than 1 /3 of the ^0.3 Mq stars that typically comprise the lens population have Jupiter-mass companions 
with semi-major axes in the range of 1.5 AU < a < 4AU. Since orbital periods of planets at these radii 
are 3-15 years, the outer portion of this region is currently difficult to probe with any other technique. 

Subject headings: gravitational lensing - planetary systems ~ stars: late type, low-mass — techniques: 
photometric 



1. INTRODUCTION 

Searches for extrasolar planets are being carried out us- 
ing several methods. More than 60 planets have been dis- 
covered by the doppler-shift technique (Marcy, Cochran, 
& Mayor 2000). While ground-based astrometric searches 
have not yielded any definitive detections, future astromet- 
ric satellites are expected to radically improve the sensi- 
tivity of this technique (Lattanzi et al. 2000). The oc- 
cultation method has yielded an important null result for 
planets in 47 Tuc (Gilliland et al. 2000), and one confir- 
mation (Charbonneau et al. 2000; Henry et al. 2000). All 
of these methods are fundamentally restricted to planets 
with orbital times shorter than the experiment, and hence 
to relatively close (and also generally massive) compan- 
ions. 

Microlensing provides a method to search for planets 
(Mao & Paczyhski 1991) that does not suffer from this lim- 
itation. When two stars are approximately aligned with 
an observer, the nearer star (the "lens") splits the light 
from the more distant star (the "source") into two im- 
ages whose brightnesses change as the relative alignment 
changes. The characteristic angular scale, 6'e ("Einstein 
ring"), and timescale, <Ej of such a microlensing event are 
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where and Ds are the lens and source distances. 
For typical events seen toward the Galactic bulge, 6'e ~ 
320^as(M/0.3Mo)i/^ which is too small to resolve di- 
rectly. However, since the images are magnified, the event 
can be identified photometrically. For typical bulge events, 
/iici ^ 25 kms~^ kpc~^, so ^e ~ 20 days, and hence nightly 
monitoring is sufficient to find most events. Four groups, 
OGLE, MACHO, EROS, and MOA have carried out such 
microlensing searches and combined have detected over 
700 events, most in the direction of the Galactic bulge 
(Udalski et al. 2000; Alcock et al. 1997; Abe et al. 1997). 
All four teams recognize these events in real time and elec- 
tronically alert the community soon after the onset of the 
event. 

If the lens has a planet that lies close to one of the lensed 
images of the source, that image is further perturbed and 
the magnification changes significantly during a time tp, 
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where M is the mass of the lens, Z^ici = AU/tTp 



where is the mass of the planet. Hence a planet betrays 
itself as a short (~ 1 day) "bump" on an otherwise nor- 
mal single-lens light curve (Refsdal 1964; Paczyhski 1986). 
Gould & Loeb (1992) showed that, with photometric pre- 
cision of about 1-2%, Jupiter-mass planets present a rea- 
sonable probability for detection throughout a wide zone 
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centered on the Einstein ring. For typical lens distances 
~ 6kpc, this sensitivity peaks at projected separations 
2 AU(Af/O.3M0)^/^. Since tp is of order or shorter than 
the sampling time of the microlensing search teams, Gould 
& Loeb (1992) advocated setting up a globe-straddling 
network of observatories to do continuous follow-up obser- 
vations of alerted events. 

The PLANET collaboration was formed in 1995 (Al- 
brow et al. 1996) expressly to carry out such observations 
and demonstrated in that pilot year that the program was 
feasible (Albrow et al. 1998). PLANET obtained substan- 
tial observing time at four southern locations (Tasmania, 
Western Australia, South Africa, and Chile) during 1995- 
99. During these five years, we monitored 50 microlens- 
ing events sufficiently well to have good to excellent sen- 
sitivity to planets, none of which displayed a clear photo- 
metric anomaly that was best explained by a planet orbit- 
ing the lens. We quantify this statement by characterizing 
the statistical sensitivity of our five-year data set to plan- 
ets. We then use these results to build mass-separation 
exclusion diagrams for the typical Galactic stars (i.e., mi- 
crolenses) that our survey probes. The basic method of 
analysis is given in Gaudi & Sackett (2000) and applied to 
event OGLE 1998-BUL-14 in Albrow et al. (2000b). The 
details of its application to the present data set are given 
by Albrow et al. (2001). 

2. PLANET FIVE-YEAR PHOTOMETRIC DATASET 

Our data were acquired over five years from six tele- 
scopes: the Canopus Im near Hobart, Tasmania, the 
Perth/Lowell 0.6m at Bickley, Australia, the Elizabeth Im 
at the South African Astronomical Observatory (SAAO) 
at Sutherland, South Africa, the ESO/Dutch 0.9m at La 
Silla, Chile, the Yale Im and the CTIO 0.9m at the Cerro 
Tololo Interamerican Observatory (CTIO) at La Serena, 
Chile. Data were collected in Cousins / and Johnson V, 
with strong emphasis on the former. 

The bulk of events analyzed here have median sam- 
pling times of 1 to 2 hours, or O{10~^t^). Observations 
from four observatories permit round-the-clock monitor- 
ing of events; the best events rarely have gaps in the data 
longer than a day. For such events, the typical photomet- 
ric precision (as judged by the scatter) is of order 1-2% for 
points near the peak, which contain most of the sensitiv- 
ity to planets. The dataset contains 5-6 high-magnification 
events with well-sampled peaks. These events contribute 
at least 1/2 of our overall sensitivity. For the final event 
sample, the median number of photometric points within 
of the peak magnification is about 140, with 75% of 
events having more than 65 points and 25% having more 
than 250 points within of the peak. Details are given 
in Albrow et al. (2001). 

3. EVENT SELECTION 

We begin with the complete sample of Galactic bulge 
events monitored by PLANET during 1995-1999, discard- 
ing those of extremely poor quality and those known to 
contain anomalies characteristic of roughly equal-mass bi- 
naries or other anomalies unrelated to binarity. Normal 
point-source/point-lens (PSPL) microlensing events are 
described by 

Fit) ^ FsAit) + Fb, 



A[uit)] 



u(t) 




4)1/2' 



[t - to)' 



(3) 



where F is the observed flux, Fg is the source flux, Fi, is 
the flux from any background light that is not magnified, 
uq is the projected source-lens impact parameter, and to 
is the time of closest approach. Separate Fs and F^ are 
required for each observatory and wave band. 

Non-planetary light curve anomalies are identified 
through explicit modeling. We discard such anomalous 
light curves from our analysis because we do not currently 
have the ability to systematically search for planetary sig- 
natures in the presence of these secondary effects. We 
note, however, that none of the excluded light curves show 
signs of short duration bumps, except MACHO 97-BLG- 
41 for which the deviation is explained naturally by binary 
rotation (Albrow et al. 2000a) rather than a planet orbit- 
ing a binary (Bennett et al. 1999). 

To eliminate events of particularly poor quality for 
planet detection, we introduce event selection criteria: 

1) All data must pass certain quality tests: only DoPhot 
(Schechter, Mateo, & Saha 1993) types 11 and 13 are ac- 
cepted. 

2) There must be at least 20 data points from at least one 
observatory in one band. 

3) Each observatory-band dataset included must contain 
at least 10 data points. 

4) The error in uq from a combined fit must be less than 
50%. 

If uq is not well constrained, then the source's path 
through the Einstein ring is not well determined, and 
hence it is difficult to estimate the sensitivity of the event 
to planets. When available, we use OGLE and MACHO 
data to help constrain mq, but not to search for planets. 
Our final sample includes a total of 43 events. A full list of 
these events along with the photometric data is presented 
in Albrow et al. (2001). 

4. SYSTEMATIC SEARCH FOR PLANETARY SIGNATURES 

The method of measuring the sensitivity of an event to 
the presence of planets, and at the same time searching 
for planetary signatures if they are present, is thoroughly 
described in Albrow et al. (2001). Very briefly, we obtain 
a PSPL fit using equation (|^), and simultaneously renor- 
malize the photometric errors at each observatory so that 
the total XpsPL this fit is equal to the number of degrees 
of freedom. Note that outliers are not included for error 
renormalization, but are included when searching for plan- 
ets. We also include in all model fits a term that accounts 
for the correlation of the photometry with the seeing that 
we observe in most of our light curves (both microlensed 
and constant stars). Although these procedures do bias us 
against binaries, the bias is only serious if all the points 
from one observatory are concentrated in a short span of 
the light curve and there are no contemporaneous data 
from other observatories. From direct inspection of the 43 
lightcurves in our event sample we find that such bunching 
affects 0(1%) of our total lightcurve coverage, which leads 
to an overestimate of our detection efficiency of a simi- 
lar magnitude. Since this bias is an order of magnitude 
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smaller than our statistical errors, we ignore it. On the 
other hand, Monte Carlo experiments with constant stars 
reveal that error renormalization and removal of system- 
atic effects are essential in order to draw reliable inferences 
from the light curves and to avoid spurious detections. See 
Albrow et al. (2001) for details and a thorough discussion. 

For each planet-star mass ratio q and each planet-star 
projected separation 9-Ed, as well as each angle a of the 
source trajectory relative to the binary axis, we find the 
best fit to the remaining parameters {to^t^.UQ, Fs, Ft). 
The corresponding thus yields Ax^ = ^ XpsPL' 
for which we set a threshold value Ax^in = 60. If 
Ax^ > Ax^. , then the geometry {d,q,a) is excluded. If 
Ax^ < — Ax^in, we tentatively conclude that we have de- 
tected a planet. The detection efficiency ei{d,q) for plan- 
ets with this (d, q) in event i is then just the fraction of all 
angles a (out of 2-k) that are excluded. We define a "plan- 
etary system" as a binary lens with mass ratio q < 0.01. 
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FIG. 1 Exclusion diagram for pairs of planet parameters {d,q), 
where d is the projected separation in units of the Einstein ring, and 
q is the planet-star mass ratio. The inner contour indicates the (d, q) 
for which the fraction of lenses with a planet is /< 1/4 at 95% con- 
fidence. Other contours are for / < 1/3, 1/2, 2/3, and 3/4. A mass 
ratio q = 0.001 corresponds approximately to rup = 0.3-Mjup, and 
a projected separation of d = 1 corresponds to a physical projected 
separation Tp = 2AU. We have assumed a detection threshold of 
Ay^ . = 60. 

^min 

We set the threshold Ax^nin = 60 by first noting the 
continuous distribution of Ax^ <J 50 in our data. If a 
significant fraction of these were due to planets, the distri- 
bution would extend to much more extreme values, since a 
small random change in the impact parameter could easily 
increase Ax^ to several hundred. Hence, the great major- 
ity of these deviations must be due to small unrecognized 
systematic errors. Monte Carlo tests performed with con- 
stant stars reveals that deviations of Ax^ <, 60 are easily 
explained by systematic and statistical noise. We there- 
fore set the threshold high enough to exclude these non- 



planetary sources of noise. We also show results for the 
more conservative threshold Ax^in = 100. 
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FIG. 2 Exclusion diagram for planets anywhere in a continu- 
ous range of semi-major axes centered on the Einstein ring. Bold 
curves show the excluded fraction (at 95% confidence) anywhere in 
the range 1.5 AU < a < 4AU, while solid curves show the frac- 
tion for lAU < a < 7AU. We have assumed a detection threshold 
of Axjj^ijj = 60, and that the primary has an Einstein ring radius 
r'E = 2 AU. For other Einstein ring radii, the separations scale as 
(rE/2 AU). The dashed curves show the negligible effect of including 
the finite size of the source in the modeling. The dotted lines are for 



100. 



This search procedure identifies two possible candidates, 
but these are also representatives of two classes of phe- 
nomena that must be excluded from our search: nearly 
equal-mass binary lenses and global-asymmetry anomalies. 
MACHO 99-BLG-18 has a ~ 15 day anomaly of amplitude 
~ 2%. Such an anomaly is longer than that expected from 
planets with q <^ 0.01, and we therefore systematically ex- 
plored binary-lens fits with q > 0.01. This uncovered a 
fit with q 0.2 that is favored over the best-fit planet 
{q < 0.01) by Ax^ = 22. We therefore exclude MA- 
CHO 99-BLG-18 from the analysis. Although Ax^ = 22 
is below our normal threshold (Ax^ = 60), we estimate 
that the probability that we have inadvertently thrown 
out a real planetary detection that is <^ 10%, smaller than 
the statistical errors on our resultant limit on planetary 
companions. See Albrow et al. (2001) for a discussion. 
OGLE 99-BUL-36 displays an overall asymmetry that is 
consistent with a distortion caused by a g ~ 0.003 planet. 
Such parallax asymmetries are a general feature of paral- 
lax which must be present at some level in all microlensing 
events (Gould, Miralde-Escude & Bahcall 1994). Indeed, 
this distortion is equally well-fit by such a parallax asym- 
metry model. Given that OGLE 99-BUL-36 has a rela- 
tively short timescale (^e — 30 days), one might naively 
expect the parallax interpretation to be unlikely. However, 
the magnitude of the parallax asymmetry is quite small. 
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and only detectable due to the high quality of the data. 
The resulting asymmetry implies reasonable values for the 
most probable mass and distance to the lens. We conclude 
that we cannot reliably detect planets from global asym- 
metries, and should exclude from our analysis all events 
that display such anomalies and reduce our efficiency esti- 
mates accordingly. Although we do not explicitly do this, 
this results in detection cfRciencies that are overestimated 
by a negligible amount, since a very special planetary ge- 
ometry is required to produce a global asymmetry (rather 
than a short duration anomaly). This conclusion is borne 
out by explicit simulations (Albrow et al. 2001). Thus 
there are no viable planet candidates out of our original 
sample of 43 events. 

5. EXCLUSION DIAGRAM FOR GALACTIC PLANETS 

From the above analysis of each event i, we obtain an ef- 
ficiency ei{d, q) as a function of planetary geometry (d, q). 
Let f{d,q) be the fraction of lenses having a planet at 
(d, g). Then, from binomial statistics, the probability of 
observing no planets is P = 1 — HJl ~ f{d,q)ei{d,q)]. 
Note that in the limit /e <C 1 (approximately valid in 
the present study), this reduces to the Poisson formula 
P = [1 _ exp(-7V)], where N{d,q) = /(d, g) ?) 
is the expected number of detections. Thus, to a good 
approximation, fractions f{d,q) > 3/ ^^ei{d,q) can be 
rejected at 95% confidence. 

Based on this analysis, we build an exclusion diagram 
(Fig. 1) based on the sample of 42 events (excluding 
MACHO-99-BLG-18) for planet parameters (d, q). To con- 
vert d and q into physical parameters of planet mass 
and projected physical separation rp, we must estimate 
the typical mass M and physical Einstein radius Dl9b 
for the events in our sample. The majority of detected 
microlensing events are almost certainly bulge stars lens- 
ing other bulge stars (Kiraga & Paczynski 1994). If the 
lenses were drawn randomly from the bulge mass func- 
tion as measured by Zoccali et al. (2000), one would ex- 
pect the typical mass to be M ~ 0.3 M© and hence the 
typical timescale to be tE ~ 20 days. However, the me- 
dian timescale for our sample is 40 days. The difference 
is probably mainly due to a bias in our selection process. 
From equation (|l|) , a bias toward large t-^ will cause biases 
toward higher M, higher tTi-oI, and lower /ij-ei- Compar- 



ing Figures la and lb from Gould (2000), we infer that 
most of the dispersion in observed timescales is due to 
/irci and TTroi, SO the bias in terms of mass is likely to be 
modest. Hence, we adopt M ~ 0.3 M© and tTj-gI ~ 40/Ltas, 
and thus 9e ~ 320 ^as. With our convention Dl = 6kpc, 
q = 0.001 ^mp = 0.3Mjup, and d = 1 rp = 2 AU. 

In Figure 2 we present upper limits for the fraction of 
lenses with planets over two ranges of semi-major axes a 
centered on the Einstein ring. To convert from projected 
separation (Fig. 1) to semi-major axis, we integrate over 
all orientations assuming circular orbits, M = 0.3 Mq and 
DlSe ~ 2AU. We calculate efficiencies using both the 
point-source approximation and allowing for finite source 
size (Gaudi & Sackett 2000) but find that the differ- 
ence is negligible. We find that less than 1/3 of lenses 
have Jupiter-mass companions anywhere in the range of 
1.5 AU < a < 4.0 AU. These are the first significant limits 
on planetary companions of M dwarfs. 
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